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THE EVOLUTION OF A GRAVITATING, 
ROTATING, CONDENSING FLUID 

In the science of Astronomy there are many observed phe- 
nomena for which, it is frankly admitted, there are at present 
no adequate explanations, while for other phenomena the 
physical causes assigned to explain them may legitimately be 
called into question. Certain of these phenomena are inti- 
mately associated with the evolution of matter in great masses 
as in a star or in the sun. Important mechanical changes with 
which observed phenomena are intimately associated must 
continually be taking place within a star as a result of the 
enormous amount of energy it is radiating into space. 

Various phases of the mechanical problem have been dis- 
cussed by mathematicians, among whom may be mentioned 
Harzer^, Wilsing^, Samson^, Wilczynski* and Emden^. Their 
conclusions differ in some important particulars and they do 
not adequately explain observed phenomena. A further con- 
sideration of the problem seems not only desirable but to be 
an imperative need of astronomy at the present time. 

Maclaurin, Jacobi, Thomson, Poincar6 and Darwin have 
discussed the problem of figures of equilibrium of rotating 
homogeneous fluids and their stabilities. They take no account 
of probable differences of density in masses of great size nor 
of differences in angular velocities as observed in the sun. 
Professor Moulton^ summarizes their work as follows: 

"In particular, considering a series of homogeneous fluid 
masses of the same density but of different rates of rotation, 
it is shown that there is a continuous series of figures of stable 
equilibrium beginning with the sphere for zero rate of rotation; 
then, with increasing rotation, passing along a line of oblate 
spheroids until a certain rate of rotation is reached; then, 
with decreasing rate of rotation but with increasing moment 
of momentum, branching to a series of ellipsoids with three 
unequal axes, and continuing until a certain elongation is 
reached; and finally, at this point, branching to a series of 
so-called pear-shaped figures. It has been conjectured that 
if it were possible to follow the pear-shaped figures suflSciently 
far, it would be found that they would eventually reach a point 
where they would sepay^: ^. iptp^j^ffiQ^ stinct masses. 
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"From this line of reasoning it has been regarded as probable 
that celestial masses, through loss of heat and consequent 
contraction, do break up in this way often enough to make 
the process an important one in cosmogony." 

It is a conmionly accepted theory that the sun and planets 
once existed as nebulae and that they, as well as the stars, 
developed from matter in a state of wide distribution. The 
problem presented by them in this condition is that of 
separate discrete particles possessing individual velocities and 
subject to accelerations produced by the force of gravitation 
due to their mutual attractions. It is essentially the prob- 
lem of the orbital motion of n bodies, and as such must con- 
form to the laws expressed by the general equations deduced 
for such motion. 

I shall therefore regard a star, or the sun, as presenting 
the n body problem in which motions are modified and energy 
produced by the impacts to which the n bodies are subject. 

In a general discussion of the problem, it would be necessary 
to consider masses of all sizes, such as might be encountered 
in nebulae, but when the configuration of particles has attained 
to the state of a star or the sun the n bodies with which we 
have to deal, will be the smallest divisions of matter, — ^namely, 
molecules. Whatever the masses composing the sun may have 
been when it was in a state of wide distribution, continued 
impact would tend to divide them into the smallest divisions of 
matter which can exist by themselves. 

The motion of the center of mass may be considered equal 
to zero since in the present problem we are concerned only 
with the evolution of n bodies with respect to it. We are thus 
led to an immediate consideration of the three integrals of area 
which are expressed by the following equations'^: 



S/ dyi dxi\ 



(1) 
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An areal velocity is essentially a vector quantity and there- 
fore not only its magnitude but also its direction must be 
specified. An arrow in the axis of rotation of length equal to 
the magnitude of the areal velocity is its conventional vector 
representation, and its origin is usually in the plane of rotation 
to which it is always perpendicular. Its direction of pointing 
indicates the direction of rotation, since it is drawn with refer- 
ence to this rotation in such a way that the arrow points in 
the direction a right-handed screw would advance if turned 
in the direction of rotation. 

Let the areal velocity of particle i be represented by a 
vector Ai, its mass by Wi, and the product utiAi by the vector 
Q which we shall call its mass rotation. 

A vector representation of mass rotations in the problem 
of n bodies would consist of n arrows of lengths equal to the 
mass rotations of the respective particles, having the center 
of mass as common origin, and pointing outward in directions 
corresponding to the rotations which they were drawn to 
represent. Let each of the n vectors thus drawn be resolved 
into components and their sums be taken along the three coordi- 
nate axes, those components which point in the positive direc- 
tion of the axes being considered positive; those pointing in 
the opposite direction negative. When this is done the geo- 
metrical significance of the constants Cx, Cy and Cg in equation 
(1) is at once apparent. They are the respective sums of 
these components along the respective axes. Thus C^ ex- 
panded takes the form 

n 

C, = C,i + C^ + (— C,8) + (-CJ+ + C^ =XC^ (2) 

i=i 
with similar expressions for C^ and C,. 

Take the vector simi of all the vectors representing mass 
rotations and we determine in direction and amount the excess 
of rotation of the configuration and the plane in which it takes 
place. This condition is conveniently expressed by the vector 
equation®, 

Xfniri-U^ = XfniA,=X.i^ (3) 
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The vectors Cx, Cy and C^ of equation (1) are the com- 
ponents of the resultant vector represented by summation Q 
of this equation. The axes of reference may be rotated so 
that one of them coincides with this resultant. Its projec- 
tions on the other two axes will then be zero. Let the Z 
axis be brought into coincidence with the resultant vector. 
It does not follow from this that there are no projections of 
vectors along the X and Y axes, representing mass rotations. 
We interpret its meaning from equation (2) by saying that the 
sum of the positive components is respectively equal to the 
sum of negative components along these axes. 

A vector representation of the conditions expressed by 

equation (1) is given in Fig. 1, the Z axis being chosen parallel 

to the resultant vector. The vector Q', in this figure represents 

2. merely the residual of mass 

rotation of the system. The 

mass rotations in the ZX and 

YZ planes and also in the 

XY plane, neglecting C/, 

Ay may have any assigned val- 



K 




^ 



An ^ y ties whatever as represented 

by Ay, Ax and A^ respect- 
ively; but in this case there 
^^* are mass rotations in the 

opposite directions of equal 
X amounts represented by the 

vectors Ay\ Kx' and A^'. In equation (2), A^. is the sum of the 
positive terms of the second member and A^.' the sum of the 
negative terms. Their sum in the case here considered is equal 
to zero. 

In polar coordinates we may write equation (3) in the form, 



mifi 



dt 



= w,A,.= C, (4) 



for particle i. The vector A»- will be constant as long as the 
resultant force acting upon it passes through the center of 
mass^, but it will vary with its deviation. As in equation (3) 
the sunmiation may be taken for all of the particles. 



Digitized by 



Google 



Evolution of a Gravitating, Rotating, Condensing Fluid 9 

Let us first consider motions of rotation only, in which all 
the particles are assumed to be describing circular orbits. The 
kinetic energy of particle i may then be expressed by the 
equation, 

K.E..=i».r.'(f)'=ii^. (5) 

Hence, the kinetic energy of rotation of a particle varies 
directly as the square of its areal velocity and inversely as the 
square of its distance from the center of mass. Kinetic energy 
may be regarded as power to do work and, if we assume A^- to 
be constant, its power to do work increases as the particle 
approaches the center of mass. 

We may put equation (5) in the following form for unit mass 

and define the kinetic potential of rotation of a particle of unit 
mass at the point r as the work which one half its centrifugal 
force would do in moving from the center of mass to that 
point. 

Taking the derivative of Vf in the direction of r we have, 

which is the outward acting force along r due to the energy of 
rotation possessed by the particle^^. It is to be noted in par- 
ticular that the kinetic potential of a particle is intrinsic poten- 
tial and is independent of the motions of other particles in the 
configuration. Each particle may thus be considered as 
possessing a specific charge of kinetic energy of rotation analo- 
gous to the electric charge of a particle in an electrostatic 
field. 

The total amount of kinetic energy of rotation possessed 
by the system is given by the equation, 



(7) 



V 



2^ K. E.^— "^ 2mf Wi „ 2» 



(8) 



which is essentially a scalar equation giving quantitatively 
the mass rotation of the system. The distinction between this 
and the vector equations (2) and (3) is to be noted. In those 
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equations negative terms cancel equal positive terms or their 
equivalents, whereas in this equation, the squares being taken, 
the total quantity of motion of rotation is included in the 
summation. 

The path pursued by a particle in the configuration may 
not be perpendicular to the radius vector, in which case the 
kinetic potential of rotation will not be the total kinetic poten- 
tial of the particle. Suppose the actual path in the time dt 

ds 
to be ds. Its velocity will then be — , and its kinetic energy 

dt 

(ds\ ^ \ ' 

— j . Resolving this velocity into components along, 

and perpendicular to r, we have the total kinetic potential, 

v/=i-[(r-.'(f)1 

It follows from this that the particle has a kinetic potential 
of translation along r as well as a kinetic potential of rotation. 

Taking the derivative of the kinetic potential in the direc- 
tion of r we obtain the force acting outward along r, 

8V/ (c^n aA 

Evaluating A in terms of r and 0, the quantity in paren- 
thesis becomes the well-known expression for the acceleration 
along a radius vector in the case of curvilinear motion. 

Taking the summation for all the particles, we find the 
total kinetic energy to be, 

The particles under consideration exist in a gravitational 
field of force, all lines of which terminate at the center of 
mass. We may therefore define gravitational potential at any 
point in the configuration as the work which the gravitational 
force acting on a particle of unit mass would do in moving from 
that point to the center of mass, supposing the field of force not 
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changed by this transfer. The gravitational potential of the 
particle according to this definition increases with its distance 
from the center of mass. Its power of doing work in the gravi- 
tational field, as in the case of falling bodies, is increased 
when we increase the distance through which it may fall* 
For the particle i we have as the expression of its gravitational 
potential. 

It is found convenient to define the two potentials with 
reference to the center of mass since the forces tmder con- 
sideration are ftmctions of distances from this point, and the 
reaction of the forces of the respective fields transfers matter 
away from or toward it. 

We speak of the two fields in terms of potential but they 
can as well and as correctly be specified in terms of energy; 
and, in the discussion following, the terms kinetic energy and 
gravitational energy will be freely used. The particle, accord- 
ing to these definitions, may be regarded as subject to two 
fields of force, in one of which the forces act outward from the 

center of mass 1 — — ; while in the other their direction is toward 
L 8r J 

the center of mass + — I . 
L 8rJ 

Hence, if a particle have a given kinetic energy of rotation 
and of translation along r [Ar positive] ^ and be subject at the 
same time to a gravitational energy, the total energy of the 
particle will be less than if subject to kinetic energy alone. 
The force due to kinetic energy will be diminished by that of 
gravitational energy, since they act in opposite directions, and 
this diminished force acting on a given mass through a given 
distance will do less work than the greater force due alone to 
kinetic energy. We have therefore for the energy of a particle 
subject to both fields, 

V, = Vj,E,-Voi. (13) 

According to this conception of the problem, V< represents the 
difference in amount of the work which the two fields would 



Digitized by 



Google 



12 University of Cincinnati Studies 

do if they acted separately on equal masses through equal 
distances. 

Particles of the configuration for which V< is positive will 
have an excess of kinetic energy and, retaining this excess, will 
describe hyperbolic orbits with respect to the center of mass. 
Particles for which V< is zero will have equal quantities of ki- 
netic and gravitational energy and their orbits will be parabolas. 
The particles for which V< is negative will have an excess of 
gravitational energy and their paths will be ellipses. 

As a simple illustration of this conception of the problem, 
take the following simple equation employed to determine the 
orbit in a two body problem. 

ic=iy +ir^y-— • (14) 

To simplify further the illustration, suppose the motion one 
merely of translation. 
We have then, 

, (dry K^M 

in which |C==V<, i(~ ) =^keu and = Vg». 

(J \ 2 

will be greater than the gravitational energy given by the term 

; if C be zero, kinetic energy will be equal to gravitational 

r 

energy; and if C be negative, gravitational energy will be in ex- 
cess of kinetic energy. Hence, a comet approaching the sun 
describes an hyperbola, a parabola or an ellipse if its kinetic 
energy is respectively in excess of, equal to, or less than its 
gravitational energy. 

Our present problem is primarily concerned with those par- 
ticles which remain a part of the configuration and hence with 
those whose gravitational energy is in excess of kinetic energy; 
but it is of importance to note that a particle having an excess 
of gravitational energy may attain parabolic or hyperbolic 
velocity either by decreasing its gravitational energy or by 
increasing its kinetic energy. 
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Taking the summation of equation (13) for the n particles, 
we have as an expression of the total energy of the configura- 
tion, 

ft ft ft 

Xyi = XyKEi- Xvoi (16) 

which is the Vis Viva or energy integral in the problem of n 
bodies. 

Put the energy equation (13) of the particle i in the form 

Suppose there are p particles in the configuration for which 
Vi has the same value. A solution of this equation will give 
us the locus of the points occupied by them, since V< is ex- 
pressed as a point function of space coordinates. 

The evaluation of the last term in the parenthesis is a 
problem of some difficulty as it involves a knowledge of the 
distribution of the particles as effected by their kinetic en- 
ergy. The discussions of figures of equilibrium, to which 
reference has already been made, are primarily concerned with 
the solution of equation (17) and the evaluation of its last 
term in the special case of rotations with equal angular 
velocities about a common axis. It is evident, however, that 
a general solution requires us to consider rotations in all 
planes as well as possible translations along all radii. We 
may, however, discuss a special case, assuming the summation 
known, and draw certain conclusions which will apply to other 
configurations. 

Let it be assumed, therefore, that a given configuration is 
spherical and that the quantity of matter per unit volume is 
the same throughout; or, to make the case more general, that 
it increases uniformly toward the center of mass. Let us 
also impose the condition that motions are motions of rota- 
tion only, and that each particle is sustained in its orbit, 
or in its position in the configuration, by its centrifugal force. 
Particles not occupying such positions will be considered as 
subject to forces of constraint. For example, a particle 
situated upon the earth's surface is subject to such a force by 
virtue of which it is constrained to revolve about the earth's 
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axis at a distance greater than it would if subject only to its 
kinetic and gravitational energy. If the matter composing the 
earth were to offer no resistance to the motion of the particle, then 
a discussion of the energy equation, (17), would give us a knowl- 
edge of the path described and of the properties of the motion. 

If all particles in a given region be subject to the same 
force, and either by an increase or by a decrease in their re- 
spective kinetic energies this force be altered, then the devia- 
tions may be discussed or considered as forces of constraint. 

A vector representation of the mass rotation of the config- 
uration under consideration would consist of narrows with origin 
at the center of mass and pointing outward in all directions. 

Through any particle i in the configuration suppose a 
spherical surface drawn concentric with the center of mass. 
The particles effective in producing gravitational energy at 
the point occupied by it will be included within this surface. 
The particles outside constitute a spherical shell of which it 
is a property that it exerts no force of attraction upon particles 
within^^. This property of a spherical shell enables us to 
evaluate the last term in the parenthesis of equation (17) for 
every particle in the configuration under consideration. It 
may be expressed by the following equation, 

X^=f 08) 

i-1 ^^ ^^ 

in which M^ is the quantity of matter or the sum of the masses 

of the particles within the spherical surface through nt^-, and 

Ti is the distance of the particle from the center of mass. 

For a circular orbit^^, 

Hence, for the case tmder consideration, in which the velocity 
along a radius vector is assumed to be zero, the energy equation 
of the particle i takes the form 

^'=^*, or V = M,r,; (20) 

with similar expressions for other particles. 

If 'we impose the condition that particles retain their 
relative positions in the configuration under different assumed 
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conditions of rotation, the summation of this equation for the 
n particles may be put in the form 

SV^^MSn (21) 

<-l <-l 

in which M is constant and equal to the summation for M<. 

The summation for r< in the second member of this equation 
is the sum of the distances of all the particles from the center 
of mass, and it may be considered to be a quantitative expres- 
sion of the density of the configuration. 

Hence, with an increase in the areal velocities of the par- 
ticles there is a decrease in density, and with a decrease in 
their areal velocities there is an increase in density. It is 
with this latter condition that we are primarily concerned, 
since the sun and stars are supposed to be bodies or configura- 
tions of particles in a state of contraction. 

Impact of particle with particle will be one of the char- 
acteristic phenomena of the configuration, and a consideration 
of its effect upon the kinetic energy and distribution of the 
particles is a matter of importance. Suppose two particles 
of equal mass, traveling in opposite directions with the same 
velocity, to meet in direct central impact. They will come 
to rest instantaneously, and whether they remain so or resume 
their motions in opposite directions will depend upon their 
elasticity or their power of restoring the motion lost by virtue 
of intrinsic forces of reaction. When the motion is wholly 
restored we say that the coefficient of restitution of the particles 
is unity; when it is wholly lost the coefficient is zero; and 
betwfen these values the coefficient of restitution is a f imction 
of the quantity of motion lost. 

Let Vi and V2 equal the linear velocities of the particles 
before impact, v and v^ their velocities after impact; let e be 
the coefficient of restitution. The kinetic energy lost by the 
impact of two particles of masses mi and W2 is given by the 
following equation^^, 

i (mi ^ + m2 v'^) = i (mi Vi^ + ma ^2^) 

*Wi m2 2 ov 
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It will be suflScient for our present purpose to treat the 
special case in which we put 

mi=W2, Vi=i^ and v=v\ 
Equation (22) for a single particle then becomes, 

|m^=^mt;i2— ^mvi^a— e^),^ (23) 

in which the last term is the energy lost by a single particle 
through impact. 

In Fig. 2, let Wi and m2 be two particles having mass rota- 
tions represented by the vectors Ai and A2 respectively. As 
represented in the figure, they are moving in opposite directions 

in the plane XY and are sus- 
tained in their orbital motions 
about at the distance r by 
the centrifugal force of their 
2 rotations. Approaching each 
"y^ other in opposite directions as 
indicated in the figure, they 
undergo direct central impact. 

^/^ /^ As a result the vectors Ai and 

Fig. 2 A2 will decrease in magnitude 

to the value zero. If the particles possess the property of 
elasticity these vectors will then increase in value but will point 
in directions opposite to those which they had before impact; 
that is, Ai and A2 will exchange places in the figure, and their 
magnitudes will depend upon the coefficient of restitution of 
the particles. 

Consider their respective masses and linear velocities to be 
equal. Let their linear velocities before and after impact be 

A A * 

written in the forms — and — , respectively, as in equation (5). 
r r 

Substitute these values in equation (23) and we have, 

^m^^im^-^m^d-^). (24) 

and 

ft n n 

V=Vr-Ai1(l-*^).or S V=SAi1-SAi!(l-^.?). (25) 

i-1 i-1 i-1 

When e is zero, the kinetic energy of rotation is wholly 
lost. The vectors Ai and A2 (Fig. 2) will be reduced to the 
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value zero and the motions of rotation will cease. No longer 
sustained in the configuration by centrifugal force the particles, 
if subject to no constraints, will move along gravitational 
lines of force toward the center of mass. When e is unity no 
energy is lost and A after impact will equal Ai (equation 25). 
The particles will then retrace their circular orbits at the dis- 
tance r from the center of mass until impact again occurs. 

Consider a thin spherical shell of radius r made up of 
particles having areal velocities equal to Ai. Impacts will, in 
this case, be so frequent that particles will describe oscilla- 
tions instead of revolutions. The summations of equations (24) 
and (25) for an interval of time will give us a knowledge of the 
rate at which the particles are losing their areal velocities and 
the surface in consequence contracting. And in the configura- 
tion under consideration their sunmiations for all the component 
shells will give us a knowledge of the rate at which it is losing 
energy of rotation and in consequence contracting. 

The effect of impact upon the vector representation of the 
total mass rotation, as given in Fig. 1, will thus be a gradual 
decrease in the magnitudes of the vectors A,., Ay and A, with 
a corresponding decrease in the vectors KJ, Ay' and A,'. If 
these vectors eventually attain the value zero, C/ will be the 
residuum of rotation possessed by the system. Repeated 
impact, therefore, reduces the configuration to one in which 
rotations take place in the same direction about a common 
axis. 

It is known that a large percentage of the energy of motion 
lost through impact is transformed into radiant energy. It 
follows, therefore, that much of the energy radiated from the 
sun and stars is that transformed by the impacts of their com- 
ponent particles through which motion of rotation is lost to 
the system. 

CoeflScients of restitution of the particles will not all be 
the same throughout, and it is of importance to know in what 
way they are effected by position in the configuration. Where 
coefficients of restitution are least, energy of motion is being 
most rapidly transformed into radiant energy; where they 
are greatest, motion is conserved. In general, we may say that 
coefficients of restitution will be least for those particles from 
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which energy generated by impact is most readily radiated 
and hence, from those at its surface. The outer shell of particles 
constitutes a non-conducting substance of greater or less 
efficiency which offers an obstruction to the passage of radiant 
energy. 

Coefficients of restitution may, then, be said to increase 
from surface to center of the mass. The condition is exempli- 
fied by steam in a boiler under pressure. The energy of the 
steam will be radiated less rapidly if the boiler be covered 
with a non-conductor, or, in other words, the coefficients of 
restitution of the steam particles will be increased by its 
addition. The kinetic energy of their motions will be trans- 
formed less rapidly into radiant energy when an obstruction is 
offered to its free passage. 

Consider the configuration to be divided up into concen- 
tric shells by concentric spherical surfaces designating them 
As, A4, As, etc., cotmting from surface to center as in Fig. 3. 

Particles of the outer 
shell will lose their 
energy of rotation 
most rapidly because 
of the smaller values 
of their coefficients of 
restitution. Eventu- 
ally, they must at- 
tain the condition in 
which their energy of 
rotation will be less 
than that of particles 
composing the shell 
A4. If, in this condi- 
tion, they retain their 
relative positions, the 
particles of both shells 
These forces acting 




Fig. 3 



will be subject to forces of constraint, 
upon particles of the outer shell will tend to drive them into 
the inner shell, while forces of constraint acting upon particles 
of the inner shell will tend to drive them outward. 
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We may obtain an anal)rtical expression for these forces 
from equations already given. 
Put equation (20) in the form 

V=:^'_^ = 0. (26) 

rC u 

To obtain the resultant force acting upon the particle in the 

direction of r we differentiate V partially with respect to r 

and obtain, 

Jf the particle retain its areal velocity, A<, and be subject to 
the same gravitational mass, M, and at the same time be con- 
strained to describe its circular orbit at a distance + ^ from 
the center of rotation greater or less than u, the resultant 
force acting upon it will no longer be zero, as in equation (27). 
Add to Ti the term + a and we obtain as an expression of the 
force of constraint to which it is subject. 

Substitute for Mf its value from equation (27), and 

p - 20V (29) 

Let the areal velocity change. Then the particle will be 
subject to a force of constraint if it maintains its position 
in the configuration, the anal)rtical expression for which may be 
obtained from equation (26). Suppose A< to change by the 
value + B. The force of constraint is then given by, 

p^_^2B(2A.±B) (30) 

fi 
This last formula may be considered as applicable to the case 
under consideration for shell As. The total force exerted by 
the outer shell may be found by taking the summation of the 
last term for all its particles, noting that the negative sign is 
is to be used as indicating that the forces act toward the center 
of mass. Employing the positive sign in equation (29) the 
summation for particles of the inner shell will give a quanti- 
tative expression of the forces acting outward. 
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The dynamic conditions are such that eventually there will 
be a flow of particles from the inner shell through the outer, 
with a consequent subsidence of particles of the outer shell 
toward the center of mass. When this exchange has been 
accomplished, forces of constraint will disappear from the 
configuration but will reappear again after another interval of 
radiation. Observations of the surface of such a configura- 
tion of particles would show a periodic change of its surface 
phenomena in which a period of quiescence would be followed 
by a period of activity. 

The condition under which these forces of constraint may 
appear is admittedly one of forced equilibrium, and their de- 
velopment will depend in large measure upon the physical 
condition of the outer shell. They may be expected to appear 
when the vibrating particles composing the outer shell lose 
suflficient kinetic energy to attain the state of a fluid having 
viscosity and tenacity. Such a fluid would prevent the inter- 
change of particles between the shells until forces of constraint 
of magnitude sufiicient to rupture the outer shell developed. 
Quantitatively these forces of constraint will be a function of 
the resistance offered by the outer shell to forces which tend 
to rupture it. The phenomena attendant upon the exchange 
of particles between the two shells will consequently vary in 
different configurations from those in which they are hardly 
perceptible to those in which important physical changes take 
place. 

Let us examine a configuration of particles in which the 
forces of constraint have attained their maximum value and 
in which an exchange of particles between the outer and inner 
shells is about to take place. A rupture in the outer shell with 
radial motions of matter both inward and outward will be 
characteristic of the surface phenomena by which these forces 
at any given point are reduced to zero. 

With a decrease in the forces acting outward on the outer 
shell's inner surface, a wave motion will be produced in it 
which will travel through the shell as a condensation. As a 
result of this wave motion the pressure upon the inner 
shell will increase until the wave produced attains the surface. 
A corresponding increase in surface activity during the interval. 
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will result. When the wave arrives at the star's surface, which 
separates it from free space, reversal of phase will take place. 
At the same time a wave in direct phase — a rarefaction — 
will be started at the shell's inner surface with an accompany- 
ing decrease in surface activity. When the waves in opposite 
phase, starting simultaneously from the two surfaces of the 
shell, meet, interference in the wave motions will take place, 
accompanied by a corresponding variation in the decreasing 
surface activity. This interference will bring the particles of 
the outer shell to rest; and with the forces of constraint 
removed exchange of particles between the two shells will cease. 

A frequency curve of surface activity will hence give us a 
knowledge of the wave motion produced in the outer shell as 
the particles comppsing it gravitate toward the center of mass, 
and there will be a quantitative correspondence between sur- 
face phenomena and the forces of constraint which produce 
them. 

Let us suppose that certain observed solar phenomena 
may be correctly explained by the deductions here set forth. 
Solar eruptive prominences are, then, particles of matter 
ejected at a high velocity from the inner shell through a rupture 
in the outer shell. In a surface eruption of large area, as in sun 
spots, the spectroscope indicates that there is motion of matter 
both inward and outward along radii. We should expect this 
condition to be a necessary consequence of the interchange of 
particles between the two shells. 

The first branch of the sunspot curve, that from minimum 
to maximum, corresponds to the wave of condensation set up 
in the outer shell at the beginning of an interchange of particles. 
The first half of the second branch of the sunspot curve corres- 
ponds to the rarefaction beginning at the inner surface of the 
outer shell; and the interruption near the middle of the second 
branch indicates the interference this rarefaction encounters 
when it meets the wave in opposite phase which simultaneously 
started from the shell's outer surface. The last half of this 
branch corresponds to the return of the outer shell to a state 
of quiescence, following this interference. 

In stars where conditions are suitable for a great develop- 
ment of the forces of constraint, a variation of light will become 
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the important phenomenon attendant upon a readjustment of 
equilibrium. The striking similarity between the frequency 
curve of solar phenomena and that of the light variation in certain 
variable stars indicates the same physical origin. 

The earth is in reality a star which has passed through the 
stages of development we now observe in other stars. Because 
of its small size it was enabled to transform its kinetic and 
gravitational energy into radiant energy and to dispose of it 
by radiation in an interval short in comparison with that re- 
quired by the sun or a star. But in the stratified condition of 
its surface we find evidence of a certain order of development 
which the deductions here set forth purport, in part, to explain. 
Particles of matter periodically projected through the outer 
shell by these forces of constraint, have eventually developed 
successively into the superimposed strata which we now ob- 
serve. 

Thus far the effect of particles in rotation about the center 
of mass has alone been considered, but it is obvious that par- 
ticles may be moving in all directions. Under these condi- 
tions, tangential components may be considered as included 
in the summations and discussions of pure rotation. It remains 
therefore to consider radial components of velocities. 

We will continue our assumption of spherical distribution 
and put the energy equation (17) of particle i in the form. 

In the case under consideration V< represents the excess of 
gravitational over kinetic energy possessed by the particle. 

Its value in gravitational units may be expressed as follows: 

V. = -^. (32) 

Za 

The energy equation then becomes, 

and its summation for the n particles gives us the total energy 
of the system due to radial components of velocities. 
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From equation (23) we obtain an expression for the energy 
a paxticle loses by impact with another particle of equal mass 
moving in the opposite direction along a radius. 

dr 
The — without subscript indicates in this case the velocity 
dt 

after impact. 

Impose the condition that the particle i describes oscilla- 
tions along r between the points a and fc, the point b being at 
the greater distance from the center of mass. Suppose further 
that its reversal of direction at these points is the result of 
impacts with particles nta and mi, respectively. Imagine a 
plane drawn perpendicular to r at the point a and consider 
the forces which tend to carry particles through it with chang- 
ing conditions of impact. 

When the two particles, m^ and ffia, of equal mass, moving 
in opposite directions along r with the same velocity, impinge 
upon each other at the point a, they will rebound with no 
tendency in either to pass through the plane at a if the in- 
trinsic forces of restitution are the same in each. But let the 
forces of restitution of particle nti be decreased below those 
of fWo, then the latter, after expending force sufficient to re- 
verse the motion of ffii, will still possess kinetic energy whose 
force acting outward along r will carry it through this plane 
if its motion be unobstructed by other forces. We may 
properly designate this a kinetic force of constraint acting 
outward along r upon particle ma impinging outward upon this 
plane. 

With a diminution of intrinsic forces of restitution and 
a consequent loss of velocity and kinetic energy, the particle 
mi will become more subject to the gravitational force acting 
on it as seen from equation (33). This condition is exemplified 
in the special case in which we consider the coefficient of resti- 
tution to be zero. Under these conditions, when impact occurs 
at a, the particle is brought to rest, but with no intrinsic 
energy tending to move it away from the plane through a. It 
will be wholly subject to the gravitational field and hence 
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will be acted upon by a force which tends to carry it through 
the plane toward the center of mass. We may properly des- 
ignate this a gravitational force of constraint acting along r, 
tending to move particles through the plane toward the center 
of mass. 

We may extend the surface through a into a spherical sur- 
face concentric with the center of mass and draw the same 
conclusion with regard to all particles which undergo impact 
at points in it. Let us consider this surface to be that which 
separates the outer shell A5 of Pig. 3 from the inner shell A4. 
As already explained, coefficients of restitution are least for 
particles of the outer shell and increase toward the center of 
mass. Hence conditions are such that particles of the inner 
shell A4, having radial motions, develop kinetic forces of con- 
straint tending to move them into the outer shell, while par- 
ticles of the outer shell develop gravitational forces of constraint 
tending to drive them into the inner shell. In considering 
the phenomena of exchange of particles between the two shells, 
these forces of constraint must be added to those due alone to 
rotation which have already been considered. 

We may now properly draw some conclusions with regard to 
the source of radiant energy emanating from a star or from the 
sun. The amount of kinetic energy transformed by the par- 
ticle i into radiant energy at a given direct central impact is 
given by the term, 

from equations (24) and (34). With a loss of kinetic energy 
of rotation the particle moves toward the center of mass, and 
the gravitational force to which it is subject does work upon it 
in this motion which is eventually transformed into radiant 
energy. 

For spherical distribution we may put the last term of the 
energy equation (17) in the form, 

where Ff is the gravitational force to which the particle is 
subject. Let {r\-r^ be the distance the particle falls toward 
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the center of mass because of its loss of kinetic energy at a 
given impact. The stmimation of the term / w<F<dr< between 

the limits ri and r 2 will then give the radiant energy contributed 
by the particle in its fall. 

Let N< be the number of impacts in a unit of time. The 
quantity of radiant energy, R, contributed by all the particles 
in this interval, or the rate at which radiant energy is produced, 
is given by the equation, 

R= S f(l-^) [(f )'+ ^]^i + i^<fm ^<dr, (37) 

Not all impacts will be direct central impacts. At each 
oblique impact, however, there will be a component of force 
centrally directed, to which the deductions for the case of direct 
central impact will apply. We may therefore consider the 
STimmations of equation (37) to include the cases both of central 
and oblique impacts. 

From equation (37) it is seen that there are two processes 
by which the radiant energy of a star is maintained. One, 
given by the first term of the second member of the equation, 
is the transformation of the intrinsic kinetic energy of the 
star's constituent particles; the other, given by the second 
term, is the transformation of gravitational energy first into 
kinetic energy then, through impact, into radiant energy. 
Particles will not fall toward the center of mass unless they 
give up the kinetic energy by virtue of which they maintain 
their positions in the configuration. Hence, in order that 
radiant energy may be produced by the fall of particles toward 
the center of mass — or, in other words, by contraction — 
there must first be a transformation of intrinsic kinetic energy 
which will result in decreasing the forces acting outward 
upon the particles to maintain them in their respective posi- 
tions in the configtiration. According to Helmholtz, the 
sun's radiant energy is maintained by contraction. It is 
evident, however, from this discussion that transformed 
intrinsic kinetic energy constitutes an important part of solar 
and stellar radiant energy. 

One of the interesting and significant facts of stellar ob- 
servation is that a large number of stars change in brightness. 
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Variable stars are of several classes which differ from each 
other in length of period, quality of light emitted, exactness 
of periodicity and light variation; but they merge into each 
other with differences so slight that a sharp separation cannot 
be made between them. This fact of itself is suflScient to sug- 
gest that there are corresponding physical conditions in the 
stars which accotmt for their light variations and for the 
gradual deviation of these variations from one class to another. 
We may therefore profitably inquire into the circumstances 
which, in a group of n particles in mutual orbital motion 
undergoing impact, would produce light variation. 

The rate at which a star is producing radiant energy is 
given by equation (37). When R varies, it must be with re- 
spect to some quantity in the second member of the equation. 
Hence, a discussion of the quantities which in this equation may 
vary with R, will give us a knowledge of the physical changes 
taking place in the star. The term in this equation giving 
the amount of kinetic energy is a slowly decreasing quantity 
so that R eventually decreases to the value zero as a result 
of the loss of this energy; but it is a quantity which varies so 
slowly that it can be left out of account in a discussion of 
variable star phenomena. 

When coeflScients of restitution decrease, R increases. 
Then the equation, 

-^ = ~V?f?i (38) 

de ;fj dei 

represents the increase in the star's rate of transforming 
kinetic into radiant energy, due to a decrease in coeflScients 
of restitution. 

We have shown in a previous section that coeflScients of 
restitution are least at the surface and increase toward the 
star's center. Hence the condition under which they will 
decrease in value is that they be brought to the surface from 
interior points. When transferred in sufficient quantity they 
will produce variation of light. The type of star which depends 
for its variation of light upon the interchange of matter between 
surface and interior has already been discussed. It was shown 
what forces develop in a contracting globe to bring interior 
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particles to the surface, and we here show why this transfer 
produces light variation. It may be well to note in passing 
that the only cause of light variation in this type of variable 
is due to the variation of coeflScients of restitution. That is, 
in equation (37) €{ varies while all other terms remain constant. 
We now examine equation (37) to see what other terms in 
its second member will give a corresponding change in R, and 
it is at once apparent that the rate of producing energy will 
vary with the number of impacts in a unit of time. Thus the 
equation 

dR = ydR< (39) 

dN p^dNi 

gives us the statement of a condition under which a star may 
vary in brightness. 

Neither of the quantities — or — will be constant with 
^ de dN 

respect to the time; so in formulating a general equation to 
include the light variations of all stars it will be necessary 
to consider these terms as varying irregularly with the time. 
They must hence be evaluated for each case to be considered. 
The following equation then gives the law of light variation, 

dL 

— , of any star whose light variation is caused by a change 

in the rate, R, of producing and radiating energy. 
dL^^/dRdRY 
dt \ de dW 

The process here outlined by which particles are periodically 
ejected from within a star through its surface, yielding up their 
kinetic energy as radiant energy in the process, is one of great 
importance to astronomy, if well founded, and we shall, there- 
fore, examine in detail certain observed phenomena of light and 
velocity variation which it purports to explain. 

When an eruption occurs at a point on the surface of a 
star by which particles are ejected with velocities such as are 
indicated by the spectroscopic observations of novae and vari- 
able stars, and in quantity suflScient to produce the corre- 
sponding variations of light observed, the forces producing the 
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eruption will react in equal amount through the star's fluid 
interior upon the star's surface diametrically opposite. K the 
star is synmietrical with respect to the diameter in question, 
the physical conditions, as regards resistance to eruptive force, 
will be the same at both of its extremities. In this case 
eruptions will occur at opposite extremities of a diameter and 
particles will be ejected in opposite directions from the star in 
nearly equal quantity and at the same velocity. The light 
from such a star would emanate from three distinct sources; 
namely, from the star itself and from the two configurations 
of particles ejected in opposite directions from it. If the direc- 
tion of projection have components in line with the earth, the 
spectroscope would indicate a motion of matter both away 
from, and toward, the observer. 

Some well defined examples of this type of line shift are 
contained in the observations of novae. Because of the vio- 
lence of the eruption as indicated by the activity attained in 
the case of Nova Aurigse, and by its spectrum, the configuration 
of particles ejected toward the observer was preceded by masses 
of particles from the star's outer surface, and light from this 
source gave a dark line spectrum, while light from the other 
configuration, ejected away from the observer, coming from 
particles having the greatest account of kinetic energy gave a 
bright line spectrum. 

In studying the spectrum of a variable star, the astronomer 
chooses those lines which are sharp, well defined and of 
suflScient magnitude to produce a good image for meas- 
urement. Particles ejected in the direction of the observer will 
give light such as to produce these lines, and the spectroscopic 
observations of variable stars to determine their radial veloci- 
ties will give the velocities with which matter is ejected from 
the star toward the observer. In this case the greatest positive 
radial velocity observed will be the actual radial velocity of 
the star, and deviations from this velocity will be the velocity 
with which matter is ejected from the star toward the observer. 
The greatest positive velocity observed will occur at a time of 
minimum light and the curve indicating velocity of ejection 
from the star should be parallel to the curve of light variation. 
On the accompanying page I have rearranged the light and 
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velocity curves of V Aquilae as published by Wright from a 
discussion of Campbell's observations, to correspond with 
these deductions. [The Orbit of V Aquilae: Ap. J. Vol. 9, 
p. 59.] We take the origin of coordinates at the time of mini- 
mum light variation and greatest observed positive velocity, and 
plot the time intervals as abscissae, and the increments in light 
and decrements in positive velocity as ordinates. The curve as 
published is a smooth curve at B, being drawn through the mean 
positions of the plotted points. In rearranging the curve at 
B, I have passed it through each of the plotted points. Or- 
dinates representing the distance between the maximum and the 

rjf-axis of each curve have 
been made equal. The first 
branch of these curves,that 
from O to A, may be con- 
sidered as representing the 
condensation of the star's 
outer shell which sets in 
when eruption begins ; that 
part of the second branch 
from A to B, the rarefac- 
tion following it; and the 
interruptions at B are 
evidences of interference in the wave motion. The similarity 
of these curves to the sunspot curve is apparent and the 
parallelism which they exhibit is as evident as that of frequency 
curves of different observed solar phenomena. 

The plotted points which coincide with the velocity curve 
are no less interesting and confirmatory of these deduc- 
tions than those which deviate from it. Starting from the 
point of minimum, the velocity curve as drawn in this figure 
passes through every plotted point from m to n. From n to o 
they deviate from the curve as here drawn and no doubt are 
indicative of secondary wave motions in the star. If instead 
of drawing a smooth curve, one be passed through each plotted 
point, this section of the curve will be not dissimilar to the 
light curves of certain other variable stars. 

With increase in the tenacity of the outer shell eruptions 
will become more violent; and this condition will be indicated 
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by an increased inclination of the first branch of the light and 
velocity curves to the X axis, with more evidence in the second 
branch of interference phenomena. 

If the globe's outer surface is not a definite fluid stirface 
but merges into outer space by imperceptible gradations, 
there will be no reversal of phase in the wave motion, and the 
rarefaction, starting from the shell's center following a con- 
densation, will continue as such. The second branch of the 
light curve will then be a smooth curve, and the forces present 
tending to restore the motion in periodic time will make the 
slope of the first branch of the curve equal to that of the 
second. 

A configuration of particles may be expected to radiate the 
same quantity of energy during successive intervals of time. 
This energy represents kinetic energy lost to the system, with 
a corresponding development of forces of constraint which pro- 
duce eruptions. Periodicity in stellar phenomena is hence 
inevitable, but the physical conditions are peculiarly suited to 
bring about slight deviations from exact periodicity both in 
time and quantity. In general, we may say that deviations 
from exact periodicity will be a ftmction of the physical con- 
ditions which produce interference, increasing with these 
conditions. 

Let us postulate star systems similar to the solar system 
each composed of several stars, one or more of which are in 
periodic eruption with corresponding light variations; then by 
the principles here outlined we are able to account for all cases 
of variations in the radial velocities of stars and for light 
variation in variables. It would be diflScult to distinguish 
those cases of variation which might be due to orbital motion. 

It is evident that the energy equation of a particle i ejected 
from within a star will change greatly at the time of an eruption. 
In the energy equation, 

V. = Vkk-V« (41) 

Yi before eruption is negative; that is, the particle possesses 
gravitational energy in excess of kinetic energy. When the 
forces of constraint have developed to such an extent that an 
eruption occurs, the forces of constraint in the outer shell, 
acting through the distances its particles move toward the 
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center of mass, do work which is expended in giving kinetic 
energy to particles ejected. In terms of the equation, this work 
is expended in increasing the value of V^; and those particles for 
which its value becomes positive leave the star from which they 
were ejected with a velocity which will carry them indefinitely 
into space. Suppose k particles ejected at an eruption, for which 
the value of V< becomes positive. The problem with which we 
are then concerned becomes the fe-body problem as we trace 
the development of the k bodies into a star. 

The large radial velocities observed in solar prominences, the 
differential radial velocities in variable stars together with 
their light variations are indicative of the high velocities 
ejected particles may attain. And when we consider the 
possible conditions of surface tension which may develop, with 
forces of constraint corresponding, the conclusion seems inevi- 
table that the process is not only possible but of conmion occur- 
rence and of fundamental importance to astronomy. 

Eruptions will be confined to equatorial regions in a rota- 
ting star, and when it has attained a figure in which one equa- 
torial diameter is greater than the others, eruptions will take 
place at opposite extremities of this diameter. The statement 
of this fact is essential to a continuation of the present dis- 
cussion. Its proof will be given in a later section. 

The state of a star's rotation at the time of eruption will 

be a matter of impor- 
tance. Pig. 4 repre- 
sents a star in equa- 
torial cross section, 
in which D is the 
direction of rotation, 
V the vector velocity 
of ejection at A and 
B at the beginning of 
eruption, and v' the 
vector velocity of 
ejection at these 
points after an inter- 
val of rotation. 
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Let an eruption occur at A. The force with which matter is 
ejected at this point will react in equal amount through the 
fluid within the star upon the surface at B. Hence a feature 
of stellar eruption which attains large proportions will be 
ejection of particles in equal quantity from opposite extremities 
of a diameter and with equal velocities. Stars of equal magni- 
tude at equal distances frona the point of projection moving in 
opposite directions are hence to be expected. 

The forces of eruption which develop in a group of particles 
through loss of kinetic energy will be manifest in groups of all 
sizes. The variation of light in asteroids gives evidence of 
their presence, and activities observed in the stm, novae and 
variable stars are coordinated when explained on the basis of 
their presence. But in masses of matter equal to great num- 
bers of stars we may expect to find their activities manifested 
on a scale of corresponding magnitude. I refer to spiral and 
other nebulae. 

The nucleus of a spiral nebula is elongated. Hence, at 
successive eruptions particles are ejected from opposite ends 
of its greatest diameter; and as rotation is a condition of this 
elongation, it is evident that the direction of projection will be 
different at each succeeding eruption, and the stars formed 
from the particles ejected will arrange themselves in a two 
branch spiral in the plane of rotation. Observation is in accord 
with this deduction, for in photographs of spiral nebulae two 
branches are characteristic and **no other number has certainly 
been found^*." An examination of these photographs for stars 
of equal magnitude situated at equal distances from the nucleus 
and on a line passing through it indicates that stars thus sit- 
uated are a not uncommon feature of spiral nebulae. An 
example of such a configuration is contained in Keeler's Volume, 
Lick Observatory Publications, Vol. VIII, Plate 23. 

The significance of star streams, the deviation of their 
directions of motion from exact parallelism, and the motion of 
groups in opposite directions are at once apparent when con- 
sidered in the light of these deductions. 

The nucleus of a spiral nebula is of great dimensions and 
light variations accompanjdng its eruptions might escape ob- 
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servation; but even on this point observational data are con- 
firmatory, as variations of light have been observed in them. 

It is evident that a positive value for all particles ejected 
at an eruption will not be given to V^ of the energy equation, 
and in many cases it must remain negative for all particles. 
We are thus concerned with particles which remain a part of 
the original system after ejection and we may properly desig- 
nate the problem which they present the g-body problem. 

In order to give proper consideration to the g-body problem, 
it is essential to know in what quantity the energy imparted to 
a particle i at eruption is distributed among the terms of its 
energy equation by the circtunstances of the eruption. 

Consider particle i, Fig. 4. At the beginning of eruption 
it starts toward A, the point on the surface where the eruption 
occurs, acquiring a velocity in that direction. The particle t, 
however, because of its great distance below the surface, has a 
linear velocity in the direction of rotation less than points on 
the surface at A, so that its resultant vector velocity with refer- 
ence to A will bring it to the outer shell in a position follow- 
ing A, as at A' in the figure. Manifestations of this phenom- 
enon are evident in sunspots where a large spot appears to 
disintegrate into numerous small spots following it. But if the 
surface at A' have sufiicient tenacity it will prevent eruption, 
and the particle i will be constrained to pass through the 
opening at A, acquiring in its motion a tangential component 
of velocity. The work expended in giving tangential velocity 
to the particle reacts in such a way as to decrease the mass 
rotation of the star. 

The path of the particles is indicated in the spiral of the 
figure. The particle in this case will be ejected with the vector 
velocity v' and will have the radial and tangential components 
Vr and Vi. Hence, in the energy equation, 

the kinetic energy of rotation, the kinetic energy of translation 
and the gravitational potential will be increased. Under these 
circumstances the particle will describe an elliptic orbit about 
the star from which it is ejected, as will all particles whose 
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kinetic energy has been thus similarly increased. In sufficient 
quantity they will form a complete ring about the star from 
which they were ejected. Such a system is unstable under 
ordinary conditions, dynamically, and must, under the mutual 
gravitation of its particles, develop into a star revolving about 
the one from which it was ejected. We are thus able to account 
for the prevalence of double stars, and for the planets and their 
satellites. And if we accept the evidence of the position of 
the solar system in the disc of the milky way as indicating an 
order of development here outlined — and other observational 
data on this point are not entirely lacking — then we may say 
that the stars of the milky way together with the sun were 
once component parts of a spiral nebula from whose nucleus 
they have been projected radially in all directions in the plane 
of rotation. 

The earth and moon possess features of marked distinction, 
not the least of which is that the earth has an atmosphere of no 
inconsiderable mass while the moon has none. The observed 
mechanical difference between the two is that the earth's rate 
of rotation is rapid in comparison with that of the moon and 
we naturally look to this difference for an explanation of the 
phenomenon. 

Let us then suppose a case in which the earth be brought 
to a state of slow or zero rotation and impose the mechanical 
changes necessary for it to sustain its atmosphere in its present 
state of distribution. In order that the atmosphere be sus- 
tained by the same centrifugal force as at present, it must 
retain its present mass rotation about its axis. In order that 
this may be accomplished the mass rotation of one half its 
particles will have to be reversed in direction and so dis- 
tributed that no motion of particles en masse with respect to 
the surface occurs. Consider Fig. 1 as representing the mass ro- 
tation of the earth's atmosphere. The change in this repre- 
sentation necessitated by the condition here imposed would be 
to reverse in direction one half the vector C/. The distribution 
required would be that each alternate particle in planes 
parallel to the equator at equal distances from the axis of rotation 
retain its present state of rotation and that each other particle 
be given a rotation of equal amount in the opposite direction. 
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The circumstances of this state of rotation may best be 
understood by referring to Fig. 2. In the present state of the 
atmosphere t»i and W2 would be rotating in the same direction 
about the point O, a correct representation of which would 
require Ai and A2 to be drawn in the same direction. But in 
the condition imposed by the earth in a state of no rotation, 
they would be rotating in opposite directions as represented in 
the figure, by the vectors Ai and A2; and if in the present state 
of the earth's rotation the vector equation of its atmosphere 
be given thus, 

C/=X ^1 + S ^2> where Ai^ = A^^ (43) 

then for the case of no rotation it would be given as follows: 

C/=2;Ai+X(-A2)=0. (44) 

In bringing about the mechanical changes necessary we do 
not change in amount the kinetic energy of the earth's atmos- 
phere; but the reversal in direction of rotation of every 
other particle would increase the number of impacts to which 
particles of the earth's atmosphere are subject. From 
equation (40) it is seen that the effect of this reversal would 
be to greatly increase the value of N, so that the rate, 

— , of transforming the kinetic energy of the earth's atmos- 
dN 

phere into radiant energy would be greatly increased. The 

condition necessary would be equivalent to that imposed by 

allowing one half the earth's atmosphere to impinge repeatedly 

upon the other half with a relative velocity of about two 

thousand feet per second. Under these circumstances the 

earth's atmosphere would lose kinetic energy as radiant energy, 

and if it be not precipitated upon the planet's surface it would, 

at least, tend to that condition. 

Hence, among the conditions conducive to a star's possess- 
ing an atmosphere is that of having an excess of mass rotation 
in the same direction about a common axis; and the amount 
of its atmosphere will be in proportion to this excess of rota- 
tion about a conmion axis whereby it is able by centrifugal 
force to maintain it without excess of impact among its con- 
stituent particles. If a group of stars be arranged in the order 
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in which the mass rotation in the same direction about a 
conmion axis increases from one to the other, their spectra will 
show a corresponding increase in absorption lines. A star's 
spectrum, therefore, may be indicative of its state of rotation. 

We may now properly take up the consideration of the 
*-body problem; — i, e., the evolution of those particles pro- 
jected from a nucleus with excess of kinetic energy. From the 
nature of the eruption by which they were ejected into space 
it has been shown that they will have slightly different direc- 
tions of projection. This will have the effect of imparting to 
the configuration a mass rotation in the same direction in large 
quantity about a common axis, and hence its spectrum will be 
one exhibiting great absorption when impact between its par- 
ticles has become frequent enough to make it visible. It will 
also be moving with a great velocity by reason of its origin. 
But there will be two forces reacting to decrease its velocity: 
one the attraction of the nucleus from which it was ejected, 
the other its collision with matter in its direction of motion, 
which because of the star's great velocity may be not incon- 
siderable in quantity. 

The excess of mass rotation about a common axis will be 
reduced by impact with other bodies, and by forces of eruption 
which may develop, and this will be indicated by spectra ex- 
hibiting fewer absorption lines. We are thus able to account 
satisfactorily for the observed fact that stars exhibiting spectra 
of much absorption are moving with greatest velocity while 
those exhibiting spectra of less absorption are moving more 
slowly. 

The forces under consideration are essentially lamellar 
vector functions, since they are derivable from a potential 
function, and a full appreciation of their significance would be 
incomplete if we omitted to discuss them in relation to the 
laminae or equipotential surfaces to which they are related. 

In the energy equation (17), we need consider only compo- 
nents of rotation, and we therefore write the equation for 
particle i as follows, 



Digitized by 



Google 



38 University of Cincinnati Studies 

When the distribution of mass rotation is spherical, that is, 
when C,' of Fig. 1 is zero, the distribution of mass is spherical; 
in which case the equation, 

r = — = a constant, 
A 

is that of the equipotential surface through the particle. This 
is the vector equation of a sphere and hence equipotential sur- 
faces in a star having no rotation are concentric spherical 
surfaces as represented in Pig. 3, considering the stirfaces there 
drawn as representing circles to suit this case. 

The forces of constraint will act along radii and will develop 
in equal amotmt at equal distances from the center of mass 
along them. There will be no tendency of particles possess- 
ing great kinetic energy to concentrate at any one point or in the 
circimiference of any given sectional plane. Eruptions may be 
expected to occur at random over its stirface. There will be 
no motion of surface strata in consequence of eruption; the outer 
shell will gravitate toward the center of mass radially from all 
points as no tangential component of force will be possible in 
this distribution. The surface will therefore remain intact and, 
if it have no atmosphere, the star will retain upon its surface the 
evidence of eruptive force to which it was at one time subject. 
The circumstances of this condition of eruption are manifest 
in a study of the moon's surface. 

We have now to discuss those cases in which, in the vector 
representation of mass rotation [Fig. 1], C,' has successively 
increasing values, and we are thus led to a consideration of 
those figures of equilibrium, or equipotential surfaces, corre- 
sponding to increasing rotation about a common axis. 

These figures have been ably discussed by the authors 
cited on a previous page and I shall accept the sequence of 
figures there stated as that corresponding to the condition of 
increasing rotation. 

I wish, however, to call attention to the physical conditions 
which require the figure to pass from the oblate spheroid to the 
ellipsoid of three unequal axes and eventually to the so-called 
pear shaped figure. 

Consider a spheroid having uniform angular velocity and 
uniform density to be divided up into equipotential surfaces. 
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Pass planes perpendicular to the axis of rotation dividing the 
spheroid into laminae or planes of finite thickness parallel to 
tfie equatorial plane. 

The energy equation of a particle i in the surface of the 
spheroid and contained in one of these laminae will be of the 
form, 

V, = ^,-^-, (46) 

and so long as the figure is spheroidal the energy equation of 
all particles in the surface of the spheroid about its whole 
circumference and contained in this laminae will have the 
same value. 

To find the force acting upon each particle in a direction 
perpendicular to the axis of rotation we take the partial 
derivative with respect to r<, thus: 

or< fi Ti Ti \ u / 

Hence as the globe contracts and r^ decreases in value the force 
Ff becomes less. As F»- approaches zero, instability among the 
particles under consideration ensues. Their mutual attractions 
will draw them together along the surface upon which the force F< 
had reacted when of sufficient magnitude to hold them in their 
respective positions in the configuration. The particles will 
congregate at a point or at different points upon the equipo- 
tential surface to which they were formerly held in uniform 
distribution, thus increasing respectively their distances from 
the axis of rotation; and hence, according to equation (47), 
increasing the force F which results in the attainment of sta- 
bility of figtu-e. A summation for all particles thus similarly 
considered will indicate the increase in force attained tending 
to hold all the particles in the configuration. 

The discussions indicate that a readjustment at this point 
may result in a variety of figures, from those which are symmet- 
rical with respect to their elongated axis to those which diverge 
widely from sjnnmetry. We may then expect a divergence in 
eruptive conditions from the one considered in which equal 
stars at equal distances from the nucleus were postulated. 

All we may say in regard to the unsymmetrical figures is 
that reactions at opposite extremities of their major axis must 
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be equal. A greater mass may be set in motion at one end of 
this axis than at the other, resulting in diflEerent velocities of 
projection such as are indicated by line shifts in the spectra of 
Nova Aurigse. And in the pear shaped figure approaching 
instability, eruption would result merely in the separation of 
the two masses. An example of this extreme type of mass 
disintegration is evident in the Double Nebula H.II 28-29 
Leonis. [L. O. Publ. Vol. VIII, plate 24.] The relative masses 
and positions of the two nebulae in this photograph, the figure 
of each, with the elongation of the larger in the direction of 
the smaller, and the configuration of stars in the field are in- 
dicative of this type of origin. 

The so called figures of equilibritma may better be desig- 
nated figures of equal activity, which as a matter of fact they 
are. In the energy equation (45), whose discussion leads to 
their determination, V^ is made to equal a constant and, when 
particle i is in the surface of the configuration, the equipotential 
surface derived coincides with the surface of the fluid. But 
an equipotential surface as applied to fluids has value aside 
from that of determining their figures; it is of primary im- 
portance in giving us a knowledge of the forces — in particu- 
lar their direction — which react to give motion to their con- 
stituent particles. We may therefore properly proceed to a 
consideration of equipotential surfaces which pass respectively 
through all particles of the configuration. 

Consider a fluid of low uniform density rotating throughout 
with the same angular velocity about an axis. In Fig. 3, a 
representation of the fluid thus rotating is given in meridian 
cross section. PP' is its axis of rotation; the ellipsoid through 
Ps its botmding surface. The curves within the figure repre- 
sent concentric ellipsoids similar to that through Pg. 
Put (0 = angular velocity, 

^ = latitude, 

K^= gravitational constant, 

Oi = axial radius, 

bi = equatorial radius, 

\ = the value such that 1 + A^ = -^ , 
a = density. 
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The term X may be considered a measure of the oblateness of 
the equipotential surface so long as its figure is an ellipsoid of 
revolution. 

A discussion of the energy equation for the case here as- 
sumed leads to Maclaurin's well known formula^^ 

In this equation the term which measures the oblateness is 
independent of the mass. The homoeoid of matter between 
the surface ellipsoid and a similar ellipsoid through the particle 
i at an internal point does not create potential at that point, 
and it follows that equipotential surfaces within the fluid are 
concentric similar ellipsoids as represented in Fig. 3. 

We shall now proceed to discuss the effect of loss of kinetic 
energy and consequent contraction in altering the oblateness of 
equipotential surfaces, starting with the condition as repre- 
sented in the figure. I shall distinguish four cases: — 

1. The effect of loss of kinetic energy in changing the 
oblateness of the eqtiipotential surface through particle i of 
shell As. 

2. The effect of this loss upon the oblateness of the equi- 
potential surface through particle i of shell A4 which loses 
none of its kinetic energy. 

3. The effect of eruption in changing the oblateness of 
the equipotential surface through particle i of shell As. 

4. The effect of eruption in changing the oblateness of 
the equipotential surface through particle i of shell A4. 

I shall assume the reader to be familiar with the vector repre- 
sentation^^ employed in deriving equation (49) and discuss 
the components of force which hold the particle in its position 
of equilibrium and determine the equipotential surface which 
passes through it. 

Let i. Fig. 5, represent the particle in equilibrium upon one 
of the equipotential surfaces Si or S2. Let PP' be the axis of 
rotation. The vector F^ represents in amount and direction 
the force of gravity; F^i and F^ components, respectively, 
of F^ and F^ + F^, and F^ is a projection of F^. The vectors 
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Let the particle lose kinetic energy. 



Fai and Pa2 are components 
of centrifugal force due to 
rotation and are tangent, 
respectively, to the curves 
Si and S2. The vectors 
Ni and N2 are normals, re- 
spectively, to the equi- 
potential surfaces Si and 
S2 at the point i. 

1. Let i be the particle 
of shell As in equilibrium 
on the equipotential sur- 
face Si. The tangential 
forces in equilibrium pre- 
venting its motion along 
the curve are Fai and F^i. 
A force of constraint will 



thereby be developed which acts in the direction F^. Suppose 
Fc to be this force. Its component tangent to Si, will be added 
to F^i as represented in the figure and the particle will no longer 
be in equilibrium upon this equipotential surface. Similarly, 
all particles of shell As will develop tangential components of 
force reacting to move them toward the poles. The particle 
i will not come to rest until particles of the shell As have, by 
piling up at the poles, developed a force of reaction equal to 
this component force of constraint. 

Suppose the configuration thus required of particles of 
shell As to be brought about by the motion of particles other 
than i. We will then be able to consider the efifect at the point 
t of loss of kinetic energy in changing the oblateness of the 
equipotential surface through it. 

Assume Fai to remain constant. Then when the force of 
constraint has been added to F^, N2 will be the resultant force 
acting on the particle, and this by definition is normal to the 
equipotential surface through t. Therefore the angular dis- 
tance between the normal to the equipotential surface at the 
point i and the line of action of gravitation through this point 
decreases, and this is the condition for decreasing oblateness. 
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Therefore the oblateness of eqtiipotential surfaces of particles 
of the outer shell decreases with loss of kinetic energy. 

2. In the case of particle i of shell A4, which retains its 
kinetic energy, F^ is negative and the angle between the 
normal and the direction of gravity increases with loss of 
kinetic energy to particles of shell As. Hence the oblateness 
of the equipotential surfaces through it, and through all other 
particles of shell A4, increases as kinetic energy is lost to par- 
ticles of shell As. Hence, whatever the arrangement of equipo- 
tential surfaces in the original configuration may have been, 
they will develop into one whose equipotential surfaces increase 
in oblateness from surface to center as represented in Pig. 6. 

We may now properly proceed to an examination of a globe 
of particles with equipotential surfaces thus arranged to see 
if phenomena which this arrangement implies are in accord 
with those observed. 

The particles of shells As, A4, etc.. Fig. 3, moving under 
the action of forces of constraint, 5 F^, develop respectively 
into hemispherical segments As', A4', etc. The particles pos- 
sessing the least kinetic 
energy are 




segment As' 



contained in 
Particles in 
the segments A4', Ag', etc. 
toward the center, possess 
kinetic energy in increas- 
ing amount which is mani- 
fest at the surface of the 
configuration as increasing 
angular velocity about the 
axis of rotation from the 
poles to the equator, as 
observed in the sun, Jupi- 
ter and Saturn. 

The resistance to eruptive force increases in proportion to 
the loss of kinetic energy; hence, resistance to eruptive force 
increases from the equator to the poles, and when the figure is 
spheroidal eruptions will be confined to equatorial regions as 
observed in the sun, Jupiter and Saturn. 
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Similarly, when the ellipsoid is one of three unequal axes, 
particles losing their kinetic energy will gravitate from the 
extremities of the longer axis and the surface at opposite ex- 
tremities of the major axis will oflEer the least resistance to 
eruptive force. The points on the surface at which eruption 
occurs will thus be determined. 

We have yet to consider cases three and four, cited above, 
concerning the effect of eruption in changing the oblateness of 
equipotential surfaces. Particles expelled will possess the 
greatest kinetic energy and their expulsion from the globe's 
interior will remove the forces of constraint existing between 
consecutive shells. A discussion of the vector figure, [Fig. 5], 
for the condition as expressed by the following equation, 

will lead to the result that equipotential surfaces toward the 
center of mass decrease in oblateness at the time of eruption, 
while those toward the surface increase in oblateness. 

In particular, the equipotential surface between shells As' 
and A4' will decrease in oblateness. The force of gravity acting 
upon particle i of shell A5', will be decreased by the expulsion of 
matter from within the globe and at the same time it will 
retain its energy of rotation. A discussion of the vector figure, 
[Fig. 5], for this condition shows that particle i will be acted 
upon by a vector component of force tangent to the surface of 
separation of the two shells, which reacts to move it toward the 
equator. And similarly all particles of shell A5' will have a 
force component reacting to move them along the surface of 
separation between the two shells. This would indicate in 
the case of the earth that the land surface of the respective 
hemispheres at a former geological epoch occupied higher 
latitudes on its surface than at present. This deduction is 
supported by the evidence of glacial action to which the land's 
surface in middle latitudes was at one time subject. The 
surface of glaciated areas bears unmistakable evidence of power- 
ful horizontal components of force directed toward the equator 
to which superimposed ice strata were subject. 

Surface particles which lose their kinetic energy through 
impact and friction tend to move to that part of the surface 
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where gravitational lines of force intersect it at right angles, 
while interior particles, retaining their kinetic energy, readjust 
themselves to lines of force in the resultant field which, in con- 
sequence of changed conditions in the outer shell, are con- 
strained to deviate more and more from coincidence with 
gravitational lines. And when, through loss of kinetic energy 
to the outer shell, reversal of lines of force occurs within the 
star, eruptions are imminent. After eruption, equipotential 
surfaces within the star approach most nearly to similarity of 
figure, but their dissimilarity increases with loss of kinetic 
energy to the outer shell until eruption occurs. Coincident 
with a period of eruption there is, then, a corresponding periodic- 
ity in the relative configurations of interior equipotential 
surfaces. 
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